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(SJ ■ ABSTRACT 



Context. The main sequence binary star 61 Cyg (K5V+K7V) is our nearest stellar neighbour in the northern hemisphere. This prox- 
imity makes it a particularly well suited system for very high accuracy interferometric radius measurements. 
Aims. Our goal is to constrain the poorly known evolutionary status and age of this bright binary star. 

Methods. We obtained high accuracy interferometric observations in the infrared K' band, using the CHARA/FLUOR instrument. 
We then computed evolutionary models of 61 Cyg A & B with the CESAM2k code. As model constraints, we used a combination of 
observational parameters from classical observation methods (photometry, spectroscopy) as well as our new interferometric radii. 
Results. The measured limb darkened disk angular diameters are 0ld(A) = 1.775 ± 0.013 mas and 9ld(B) = 1.581 ± 0.022 mas, 
respectively for 61 Cyg A and B. Considering the high accuracy parallaxes available, these values translate into photospheric radii of 
R(A) = 0.665 ± 0.005 R Q and R(S) = 0.595 ± 0.008 R Q . The new radii constrain efficiently the physical parameters adopted for the 
modeling of both stars, allowing us to predict asteroseismic frequencies based on our best-fit models. 
Conclusions. The CES AM2k evolutionary models indicate an age around 6 Gyrs and are compatible with small values of the mixing 
^\ ' length parameter. The measurement of asteroseismic oscillation frequencies in 61 Cyg A&B would be of great value to improve the 

■^J" , modeling of this important fiducial stellar system, in particular to better constrain the masses. 
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1 . Introduction parallax of ths system was first measured by Bessel ( 1838}, and 

' . , _ , , , r r , it is now known with extremely high accuracy. The proper mo- 

^ ; Binary stars are numerous in the Galaxy but only few of them tion of more than 5 „ first determined b piazzi in the 

can be fully calibrated physically. They are interesting for sev- xvmth ce makes k one of the fastest moyi stars in terms 

eral reasons, among them the most important for evolutionary of displacement. Although some of this motion comes 

modeling is that the two components are coeval, and share the from ^ imit of 6 1 c to us> the ir has also a hi h radial 

g . same chemical abundances of the original material. This makes yeloci of 1Q8 ^ i ndicating that 6l Cyg is not a member of 

their modeling much easier than for single stars, as several free ^ ^ digk of Qur Qa[ The imit of 6 1 c makes it 

parameters of the models are shared between the two stars (age a northem anal of ^ numerical modeling benchmark a Cen. 

and helium content). Possibly the best example of physical call- fts , ^ ^ means ^ ^ ^ ^ eagiest low . mass 

bration of a binary system is provided by our closest neighbour, dwarf tQ resolye interferometrically . The spectral types of its 

a Cen (Thevenin et al 2003 Its two solar-type components twQ members (R5 y and K?v) ideaU lement our ious 

were analyzed in detail with the use of photometry spectrome- ofaCenA&B (G2V+ K1 V, Kervella et al. l2003il Bigot 

try, astrometry and interferometry. It resulted in excellent know 1- et aJ 2QQ6) The masses of 6 ! Cyg A & B are controversial, 

edge of the masses of the two components (within less than 1 %), whh estimates j from imatel 0J4 and .46 M 

and, through modeling of the stars, of fundamental parameters (Gorshanov et al jf^ to 0.67 and 0.59 M (Walker et al. 1995 ). 

not directly accessible to measurements, like the age or the he- wkh effectiye tem ratures of about 4400 and 4000Kj m 

hum content. Therefore, these two stars are, with Procyon, Sinus shine wift luminosities of onl 15 and 08 Lq . There is no 

and Vega, among the most important benchmark stars for new confirmed lanet around ^ dth h indications exist that 

developments of stellar structure and ^tmosphere m^els The 6{ Cyg B could host a giant planetary companion (see. Sect.O- 

Pre i™ d ^n U f " S A u (HD 201 ° 91 ' HIP U 10421 4) As dimmer versions of soJ type stars, they have magnetic^- 

and B (HD 201092, HIP 104217), the nearest stars in the north- ^ ^ tQ ^ of ^ s J\ he J g ^ feint y er 

ern hemisphere. The two stars constitute a visual binary pair with j j (Hempelmann et al. 2006 ). Their rotation periods are 

a very long orbital period (* 650 yrs), also known as Gl 820. The of ^ Qrder of 35 day ^ therefore no rotational distorsion of their 

Send offprint requests to: P. Kervella photospheres is expected. The abundances of heavy chemical el- 

Correspondence to: Pierre.Kervella@obspm.fr ements have been determined (Luck & Heiter |2005| |2006|l in 
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these stars which are found slightly metal poor (-0.2 dex), so a 
priori older than the Sun but belonging to the galactic disk. 

In Sect. [2] we detail the interferometric observations and the 
corresponding physical parameters we derive (angular diameters 
and linear radii). Together with the additional observables listed 
in Sect. [3] we propose in Sect.|4]a modeling of the two stars using 
the CESAM2k code. We finally present asteroseismic frequency 
predictions in Sect. [5] 

2. Interferometric observations 

2.1. Instrumental setup 

Our observations of 61 Cyg were undertaken in November 2006 
in the near infrared K' band (1 .9 < A < 2.3 /mi), at the CHARA 
Array (ten Brummelaar et al. 120051 1 using FLUOR, the Fiber 
Linked Unit for Optical Recombination (Coude du Foresto et 
al. 120031 . We used the F LUOR Data reduction softwar e (DRS) 
(Cou de du Foresto et al. 119971 Kervella et al. I2004al Merand 
et al. |20061 l. to extract the squared instrumental visibility of the 
interference fringes. The baseline was chosen according to the 
predicted angular sizes of 61 Cyg A & B (approximately 1.5 to 
2 mas) and the wavelength of observation, in order to obtain the 
most constraining measurements on the angular diameters. This 
led to the choice of the CHARA telescopes S 1 and E2, separated 
by a linear distance of 279 m. The calibrator stars were chosen 
in the catalogue compiled by Merand et al. (2005), using criteria 
defined by these authors (Table [TJ. They were observed imme- 
diately before or after 6 1 Cyg in order to monitor the interfero- 
metric transfer function of the instrument. 

We selected these calibrators so that their visibility measure- 
ments have comparable signal-to-noise ratios to 61 Cyg A & B. 
As they have nearly the same effective temperature and bright- 
ness as our targets, their angular diameters are also similar. It 
results that the final precision of our angular diameter measure- 
ments of 61 Cyg is limited by the uncertainties on the calibra- 
tors' diameters. We took them carefully into account in our fi- 
nal error bars, including the correlations between different cal- 
ibrated observations of the same star. The reader is referred to 
Perrin (120031 1 for a detailed description of the error propagation 
method we used. The final calibrated visibilites of 61 Cyg A and 
B are listed in Table [2] 

As a remark, using smaller (hence fainter) calibrator stars 
would have led to a larger observational uncertainty on the mea- 
sured visibility of the calibrator and would not have improved, 
or only on a marginal note, the accuracy of the angular diameters 
of61CygA&B. 

2.2. Limb darkened angular diameters 

In order to estimate the angular diameter from the measured vis- 
ibilities it is necessary to know the intensity distribution of the 
light on the stellar disk, i.e. the limb darkening (LD). We selected 
the four-parameter LD law of Claret (120001 1: 

4 

m/i(D = i-2 fl *d -j" 1 ) (!) 

The at coefficients are tabulated by this author for a wide range 
of stellar parameters (T e B, logg,...) and photometric bands (U to 
K). To read Claret's tables for the PHOENIX models, the effec- 
tive temperatures of 61 Cyg A and B were rounded to r e ff(A) = 
4400 K and T eS (B) = 4000 K, with a metallicity rounded to solar 



Table 2. Squared visibility measurements obtained for 61 Cyg A 
& B. B is the projected baseline length, and "PA" is the azimuth 
of the projected baseline (counted positively from North to East). 



MJD Star Bjm) PAT) V 2 ± cr(V 2 ) 



53540.401 


61 Cyg B 


210.98 


-5.67 


0.1714 


±0.0114 


53540.449 


61 Cyg B 


210.67 


-16.13 


0.1840 


±0.0116 


54043.169 


61 Cyg B 


204.09 


-34.24 


0.2982 


± 0.0257 


54044.159 


61 Cyg B 


205.08 


-33.05 


0.3049 


±0.0162 


54044.199 


61 Cyg B 


197.90 


-39.38 


0.2694 


± 0.0340 


54049.172 


61 Cyg B 


231.56 


52.54 


0.1386 ± 0.0092 


54049.207 


61 Cyg B 


217.61 


41.37 


0.1919 


± 0.0098 


54052.135 


61 Cyg B 


33.41 


-23.29 


0.9670 


±0.0168 


54052.171 


61 Cyg B 


32.80 


-29.16 


0.9405 


±0.0140 


54052.208 


61 Cyg B 


31.76 


-34.56 


0.9253 ±0.0144 


54055.180 


61 Cyg B 


94.08 


69.75 


0.7731 ±0.0120 


54055.222 


61 Cyg B 


82.60 


56.28 


0.8229 


±0.0191 



Table 3. Limb darkening coefficients of the four-parameter 
power law of Claret (EOOOi l for 61 Cyg A & B. 



Star 




«2 


«3 


a 4 


61 Cyg A 


0.9189 


-0.5554 


0.3920 


-0.1391 


61 Cyg B 


1.2439 


-1.5403 


1.3069 


-0.4420 



(log[M/H] = 0.0), a turbulent velocity VY = 2km/s, and a sur- 
face gravity rounded to \ogg = 5.0. From the masses we derive 
in Sect. [4] (listed in Table [5} and the measured radii (Sect. [2~4l , 
the effective gravities are logg = 4.71 and \ogg = 4.67, respec- 
tively for 61 Cyg A and B. The values of the four a, parameters 
corresponding to Claret's intensity profiles of the two stars in the 
K band are given in Table [3] 

The FLUOR instrument bandpass corresponds to the K' fil- 
ter (1.9 < A < 2.3 /mi). An effect of this relatively large spectral 
bandwidth is that several spatial frequencies are simultaneously 
observed by the interferometer. This effect is known as band- 
width smearing. It is usually negligible for V 2 > 40%, but not in 
our case as the measured visibilities are closer to the first min- 
imum of the visibility function. To account for this effect, the 
model visibility is computed for regularly spaced wavenumber 
spectral bins over the K' band, and then integrated to obtain the 
model visibility. This integral is computed numerically and gives 
the model V 2 as a function of the projected baseline B and the an- 
gular diameter #ld- A simple^ 2 minimization algorithm is then 
used to derive #ld, the 1 <x error bars and the reduced^ 2 . For 
more details about this fitting procedure, the reader is referred to 
Kervella et al. (|2003bt or Aufdenberg et al. d2003T >. 

In the infrared K band, the LD is much weaker than in the 
visible, and for instance the difference in the LD of 61 Cyg A 





61 Cyj 


r A 


Tin 


JJJtU.J / J 


61 Cyj 




o 1 n 

ZlU. 


53540.429 


61 Cyt 


;A 


210. 


53542.433 


61 Cy^ 


;A 


210. 


54043.158 


61C yi 


;A 


205. 


54043.203 


61C yi 


;A 


197. 


54044.149 


61Cyt 


;A 


206. 


54044.190 


61Cyj 


;A 


199. 


54052.124 


61C yi 


;A 


33. 


54052.159 


61C yi 


;A 


33. 


54052.193 




;A 


32. 


54052.234 


61Cyt 


;A 


30. 


54055.169 


61C yi 


;A 


96. 


54055.211 


61Cyj 


;A 


85. 



no 


5.35 


U. 1 194 ± U.UUV2 


98 


0.21 


0.1381 ± 0.0098 


89 


-11.70 


0.1322± 0.0082 


80 


-13.88 


0.1356± 0.0084 


50 


-32.51 


0.1407 ±0.0258 


66 


-39.53 


0.1882 ±0.0122 


30 


-31.35 


0.1393 ±0.0270 


81 


-38.07 


0.1792 ±0.0207 


54 


-21.29 


0.9559 ±0.0123 


04 


-27.24 


0.9239 ±0.0166 


23 


-32.48 


0.9517 ± 0.0092 


70 


-37.98 


0.9656 ± 0.0090 


60 


72.69 


0.6990 ± 0.0107 


78 


60.16 


0.7427 ±0.0105 
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Star 


m v 


m Ks 


Spect. 


8 LD (mas) 


0udk (mas) 


a(°) 


HD 196753 
HD 200451 


5.94 
7.23 


3.34 
2.90 


KOII-III 
K5III 


1.047 ±0.013 
1.485 ±0.019 


1.022 ±0.013 
1.441 ±0.019 


16.2 
12.3 



Table 1. Calibrators used for the observations. They were selected in the catalogue assembled by Meiand et al. (2005 ). The limb 
darkened (Old) and uniform disk (Sud k, for the K band) angular diameters are given in milliarcseconds (mas). The angular separation 
a between the calibrators and 61 Cyg is given in the last column, in degrees. 



and B is small. It can be quantified by the ratio B of the LD 
angular diameter #ld and the equivalent uniform disk angular 
diameter #ud- From Claret's models, we obtain /3(A) = 1.0305 
and /3(B) = 1 .0263. The difference between the two stars is 0.4% 
for Ar e ff w 400 K. To account for the uncertainty on the effective 
temperatures, we thus included a +0.2% systematic uncertainty 
on the LD angular diameters of the two stars. The LD being 
in any case small in the K band, this appears as a reasonable 
assumption. 

There is also an uncertainty on the effective wavelength of 
the instrument. As the measurements are differential in nature 
(between the scientific targets and the calibrators), this impacts 
the measurement only as a second order effect. The true effective 
wavelength of FLUOR instrument was calibrated by observing 
with high accuracy the binary star tPeg (Meiand et al. 2008 
in prep.). The relative uncertainty introduced on the angular di- 
ameters of 61 Cyg A & B is estimated to +0.1%, i.e. negligible 
compared to the statistical uncertainties (0.7% and 1.4%). Both 
the LD and wavelength uncertainties were added quadratically 
to the errors from the fitting to arrive at the final uncertainties of 
the LD angular diameters. The result of the fits are presented for 
61 Cyg A & B in Figs.[TJ The reduced^ 2 are 1.0 and 3.6, respec- 
tively for A and B. From the CHARA/FLUOR visibilities, we 
derive the following LD angular diameters: 

6» LD (61 Cyg A) = 1.775 ± 0.013 mas, (2) 

6<LD(61CygB) = 1.581 ±0.022 mas. (3) 

2.3. Angular diameter discussion 

We can compare the measured LD angular diameters to the ex- 
pected values for these stars from the surface brightness-color re- 
lations calibrated by Kervella et al. d2004bll. U sing the (B,B-L) 
relation (photometry from Ducati et al. 2002), we obtain the fol- 
lowing predicted values: (?LD(61CygA) = 1.813 ± 0.019 mas, 
and 6»LD(61CygB) = 1.704 + 0.018mas. While the agree- 
ment is good for A (within « 1 cr), the predicted size for B 
is 4 cr larger than the measured value. Both our measured an- 
gular diameter value for 6 1 Cyg B and this estimate from sur- 
face brightness are significantly different from the measure- 
ment obtained by Lane & Colavita (2003) using the Palomar 
Testbed Interferometer. They obtained an angular diameter of 
#ld(61 CygB) = 1.94 ± 0.009 mas, that corresponds to a dif- 
ference of +8 cr from our CHARA/FLUOR measurement, and 
+ 11 cr from the surface brightness predicted values. One should 
note that the measurement of the angular size of the star was 
not the primary focus of their work (in particular, no spatial fil- 
ter was used for this measurement, resulting in a more difficult 
calibration of the visibilities). 

These discrepancies combine with the larger observed dis- 
persion of the CHARA/FLUOR measurements (compared to A) 
to indicate a possible additional contributor in the interferomet- 
ric field of view around B (~1"). Moreover, 61 CygB shows a 
slight visibility deficit at our shortest baseline which corresponds 



to a 1.0 + 0.4% photometric excess. These observations can be 
explained if B is surrounded by a disk, or by the presence of a 
faint companion. Interestingly, an 8 Mj companion was proposed 
by Strand ( 119431 119571 1, based on astrometric measurements. 
This possibility was suggested again by Deich & Orlova dl9771 l, 
Deich d 19781 1, and recently by Gorshanov et al. {2006), based on 
astrometry of the two stars. One should note however that a giant 
planet with a mass around 10 My and the age of 61 Cyg will be 
extremely faint, according for instance to the models by Baraffe 
et al. (120031 1. with an absolute K band magnitude of M# * 30, 
and an apparent magnitude of iuk ~ 28. It thus appears un- 
likely that such a planet can influence the CHARA/FLUOR in- 
terferometric measurements. Moreover, using radial velocime- 
try, Walker et al. (1995) excluded the presence of low mass 
companions around both components of 61 Cyg down to a few 
Jupiter masses, up to periods of about 30 years. 

As our present data set is too limited to investigate this pos- 
sibility, we postpone this dicussion to a forthcoming paper. The 
dispersion of the V 2 measurements observed on B is in any case 
taken into account in the bootstrapped error bars of the LD angu- 
lar diameters and will therefore not affect our modeling anaysis. 

2.4. Linear photospheric radii 

We assumed the following parallax values for the two stars: 

tt(61 Cyg A) = 286.9 + 1.1 mas, (4) 

tt(6 1 Cyg B) = 285 .4 ± 0.7 mas. (5) 

They are taken respectively from van Altena et al. 1995} and 
the Hipparcos catalogue (ESA 1997). As a remark, van de 
Kamp (119731 ) obtained slightly different values (jta - 282.8 ± 
1.3mas and7TB = 288.4+1. 6mas)using photographic plates col- 
lected over the 1912-1972 period (see also van de Kamp [19531 . 
From the reprocessing of the original Hipparcos measurements, 
Van Leeuwen (2007a, 2007b) obtained the following parallaxes: 
tta — 286.8 +6.8 mas and Kb - 285.9 + 0.6 mas, compatible with 
our assumed values. 

From the combination of the limb darkened angular diame- 
ters and trigonometric parallaxes, we derive the following pho- 
tospheric linear radii: 

R(6 1 Cyg A) = 0.665 + 0.005 R Q , (6) 

R(61 CygB) = 0.595 + 0.008 R . (7) 

The relative uncertainties on the radii are therefore ±0.8% and 
+ 1.4% respectively for A and B. Thanks to the high precision of 
the parallaxes (0.38% and 0.25%), the radius accuracy is limited 
by the precision of the LD angular diameter measurements. 

3. Additional observational constraints 

3.1. Masses 

In spite of the proximity and large semi-major axis (a = 24.65") 
of 61 Cyg, its very long orbital period of 6-7 centuries makes it 
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baseline (m) baseline (m) 

Fig. 1. Visibility data and adjusted limb darkened disk model visibility curve for 61 Cyg A (left) and B (right) (see also Table|2]i. 



Table 4. Apparent magnitudes of 61 Cyg A and B from 
ESA (1997) and Ducati (2002). The uncertainties that were not 
available were set arbitrarily to +0.05 mag. 



Band 


61 Cyg A 


61 Cyg B 


U 


7.48 ± 0.06 


8.63 ± 0.05 


B 


6.27 ± 0.03 


7.36 ± 0.02 


V 


5.20 ± 0.03 


6.05 ± 0.01 


R 


4.17 + 0.06 


4.87 ± 0.05 


I 


3.52 ± 0.06 


4.05 ± 0.05 


J 


3.10 + 0.05 


3.55 ± 0.02 


H 


2.48 ± 0.08 


2.86 ± 0.04 


K 


2.35 ± 0.03 


2.71+0.05 


L 


2.28 ± 0.06 


2.61+0.05 



very difficult to determine accurate dynamical masses. Quoting 
Zagar d 19341 1, Baize ([1950) reports a period of 692 years, while 
Cester, Ferluga & Boehm ( 1988) quote a period of 720 yrs and 
equal masses around 0.5 +0.1 M for both components. Van de 
Kamp (1 1954| l proposed masses of 0.58 and 0.54 M for 6 1 Cyg A 
and B, respectively. From the orbit determined by Worley & 
Heintz dl9831 l. Walker et al. (119951 1 quote masses of: 



M(61CygA) 
M(61CygB) : 



: 0.67 M , 
0.59 M . 



(8) 
(9) 



While these values are rather uncertain, we will use them as 
our starting guesses for the CESAM2k modeling. As a com- 
parison, the masses proposed by Gorshanov et al. (2006)) of 
M(61CygA) = 0.74 + 0.13 M Q and M(61CygB) = 0.46 ± 
0.07 M have a larger ratio, although the total mass of the sys- 
tem (M A+B = 1.20 + 0.15 M Q ) is compatible with Walker et al.'s 
value. 



3.2. Effective temperature 

Knowing the LD angular diameter of the stars, it is possible to 
invert the surface brightness-color (SBC) relations calibrated by 



Kervella et al. (2004b) to retrieve the effective temperature of a 
star from any of its apparent magnitudes. Applying this method 
to 61 Cyg A and B (Fig. 13 yields consistent effective tempera- 
tures for the BVRI bands of: 



r eff (61 Cyg A) = 4400 ± 100 K, 



r eff (61CygB) = 4040 + 80 K. 



(10) 



(11) 



The apparent magnitudes of the two stars (Table @|i were taken 
from the Hipparcos catalogue (ESA 1997) and the catalogue 
compiled by Ducati (2002). Thanks to the proximity of 61 Cyg, 
we can neglect interstellar extinction. As a remark, the average 
effective temperature retrieved from the HKL bands are higher 
than from the visible bands, respectively by 80 and 200 K for 
61 Cyg A and B. It may be a consequence of the activity of the 
stars, but more likely, it could be due to the relatively poor pho- 
tometry available in these bands. The two stars are in particular 
too bright to have accurate magnitudes in the 2MASS catalogue 
(Skrutskie et al. 120061 1. For this reason, we consider the T e $ es- 
timates based on the visible photometry (BVRI) more reliable. 
The temperatures derived from the J band are however in good 
agreement with these average values. The derived temperatures 
are identical, within their statistical uncertainties, to those com- 
puted from the temperature scale of Alonso, Arribas & Martinez- 
Roger (TT996l >: 4340 K and 3980 K. 

For comparison, Luck & Heiter (120051 1 determined higher 
values of r eff = 4640 K and 4400 K for 61 Cyg A and B from 
spectroscopy. More generally, different techniques result in sig- 
nificantly different values for 61 Cyg A and B, and discrep- 
ancies as large as 300K are found in the literature. This may 
be a consequence of the activity of the two stars (see e.g. Hall, 
Lockwood & Skif f|2007l l. 61 Cyg A is also classified as a BY Dra 
type variable, showing a low amplitude rotational modulation 
of its light curve due to star spots and chromospheric activity 
(Bohm-Vitense |2007l 
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4800 , 
4700 - 
4600 - 
r 4500 - 
r 4400 - 
- 4300 - 
4200 - 
4100 - 
4000 - 



Table 5. Observational constraints (upper part), and best-fit pa- 
rameters derived from our CESAM2k modeling of 61 Cyg A & 
B (lower part). See text for the corresponding references. 



Pa r?i m pfpr 


61 Cmp A 


61 TvoR 


n (mas) 


286.9 ± 1.1 


285.4 ± 0.7 


[Fe/H] 


-0.20 ± 0.10 


-0.27 + 0.19 


R(R ) 


0.665 ± 0.005 


0.595 ± 0.008 


Teff (K) 


4400 ± 100 


4 040 ±80 


L(L e ) 


0.153+0.010 


0.085 + 0.007 


Initial He content Y„„ 


0.265 


0.265 


Initial [Z/H] (dex) 


-0.10 


-0.10 


Final [Z/H] (dex) 


-0.15 


-0.15 


Age (Gyr) 


6.0 ± 1.0 


6.0 ± 1.0 


Mass (M Q ) 


0.690 


0.605 


a (MLT convection) 


1.2 


0.8 



Fig. 2. Effective temperature estimates for 6 1 Cyg A (top) and 
B (bottom) based on their apparent magnitudes and the surface 
brightness-color relations by Kervella et al. (2004b). 



3.3. Luminosity 

The luminosity of the two stars of 61 Cyg can be computed us- 
ing two methods: (1) combining the effective temperature and 
the interferometric radius, (2) from the apparent magnitude, the 
bolometric correction, and the parallax. 

From the effective temperatures listed in Sect. 13.21 and the 
radii measured with CHARA/FLUOR (Sect. |2~4l i, a straightfor- 
ward application of Stefan-Boltzmann's law (L = AnR 2 crT^ s ) 
gives L(61CygA) = 0.149 + 0.006 L , and L(61CygB) e = 
0.085 + 0.004 L@. 

We can also use the K band bolometric corrections (BC#) 
from Houdashelt, Bell & Sweigart (2000) to retrieve the bolo- 
metric magnitude from the apparent K band magnitudes of the 
stars (Table 0). An interpolation of Houdashelt et al.'s tables 
gives BC K (A) = 2.152 and BC K (B) = 2.433. Considering the 
uncertainty on the metallicity of the two stars, and the differ- 
ence observed between the visible and infrared effective tem- 
perature estimates (Sect. 13.2b . we attributed a +0.07 mag er- 
ror bar to these bolometric corrections. Together with the par- 
allaxes, they correspond to absolute bolometric magnitudes of 
M bol (A) = 6.80 + 0.08 and M bol (B) = 7.42 + 0.09, and the fol- 
lowing luminosities (with Mbd(O) = 4.75): 

L(61 Cyg A) = 0.153 + 0.010 L , (12) 



3.4. Metallicity 

An error of 0. 10 dex on the metallicity translates into an error of 
160 K on the effective temperature, so a good knowledge of the 
metallicity is important. 

A number of metallicity estimates can be found in the lit- 
erature. Zboril & Byrne d 19981) determined the metallicity of 
61 Cyg A to be [M/H] = -0.3 ± 0.2, from high-resolution spec- 
troscopy and LTE atmosphere models. They also found a mean 
metallicity index [M/H] = -0. 2 for a sample of 18 K and M 
field stars. Tomkin & Lambert (1999) obtained relatively low 
iron abundances of [Fe/H] = -0.43 ± 0.10 for 61 Cyg A and 
[Fe/H] = -0.63 + 0.10 for 61 CygB, based on equivalent fine 
widths and LTE model atmospheres. Using a similar method, 
Afferetal. d2005l > measured [Fe/Hl = -0.37 +0.1 9 for 61 Cyg A. 

For the present study, we choose to rely on the work by Luck 
& Heiter (|2005), who determined: 

[Fe/H] (61 Cyg A) = -0.20 + 0.10 (14) 

[Fe/H] (61 Cyg B) = -0.27 + 0.19 (15) 

For the abundance analysis, they used plane-parallel, line- 
blanketed, flux-constant, LTE, MARCS model atmospheres. 
These models are a development of the programs of Gustafsson 
et al. (119751 1. We selected these values, as they are well in line 
(within 1 to 2 <x) with the measurements obtained by other au- 
thors. In addition, we used the MARCS models for the atmo- 
sphere description of our evolutionary modeling (Sect. 3J, so 
Luck & Heiter's determination appears as a natural choice for 
consistency. 



L(61 Cyg B) = 0.085 + 0.007 L . (13) 

These figures are in excellent agreement with those computed 
from the effective temperature and the radius. In our subse- 
quent modeling of the two stars, we will use simultaneously 
as constraints the different available observables. However, var- 
ious correlations exist between these observables, for instance 
between R, L and T e ff through Stefan-Boltzmann's law, and we 
prefer to limit the influence of these correlations as much as pos- 
sible. However, this was not possible with all parameters as, for 
instance, the parallax plays a role in the derivation of most of 
them. As the bolometric corrections have the advantage of being 
independent from the radius measurement, we retain them for 
our modeling of the two stars. 



4. Modeling with CESAM2k 

The observational constraints and parameters used to construct 
our CESAM2k evolutionary models (Morel 119971 Morel & 
Lebreton 2007) are summarized in Table [5] 

As a binary system, components A and B of 6 1 Cyg must 
have the same age, initial helium content, and metallicity (as- 
suming that the system formed as a binary). By comparison with 
our previous modeling efforts (see e.g. Thevenin et al. [2005), 
we implemented a few modifications for the present work. The 
equation of state now includes Coulomb corrections due to low 
masses of these stars and the corresponding high density. We 
used MARCS models (Gustafsson et al. 120031 120071 [200W1 for 

1 http://www. marcs, astro, uu. sej 
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the atmosphere description. The treatment of convection makes 
use of a new modified routine (conv_a®) based on the classical 
mixing length theory (hereafter MLT). The mixing length van- 
ishes at the limit between convective and radiative zonesB Note 
that our small values of the a parameter (that describes the con- 
vection in MLT) are in good agreement with the systematic trend 
of this parameter in low mass stars, where its value decreases 
steeply with mass. For instance, if for the Sun the most com- 
mon quoted value of a is about 1.9 (Morel et al. 1999 ), Yildiz et 
al. (|20O6 ) obtain for low mass stars values such as 1 .2 for 0.8 M 
or 1.7 for 0.9 M Q . Compared to these values, our values of a of 
0.8 for 6 1 Cyg B (0.6 M ) and 1 .2 for 6 1 Cyg A (0.7 M G ) appear 
reasonable. 

We define the stellar radius of a model as the bolometric one, 
which is equivalent to the interferometric definition of the limb 
darkened angular diameter. The microscopic diffusion of chem- 
ical species is taken into account according to Burgers (1969), 
using the resistance coefficients of Paquette et al. £[986). 
According to the prescription of Morel & Thevenin (2002), we 
introduce an additional mixing parameterized by Rev « 1 . This 
parameterization is not important for the stars 61 Cyg A & B 
because the diffusion is inefficient for low mass stars with solar 
metallicity. 

The adopted metallicity [Z/X], which is an input parame- 
ter for the evolutionary computations, is given by the iron abun- 
dance measured in the atmosphere with the help of the follow- 
ing approximation: log(Z/X) ss [Fe/H] + log(Z/X) . We use 
the solar mixture of Grevesse & Noels ( fT993b : (Z/X) Q = 0.0245. 
The evolutionary tracks are initialized at the Pre-Main Sequence 
stage ; therefore, the ages are counted from the ZAMS. 

To fit the observational constraints (r e ff , L and surface metal- 
licity [Z/X] sur f) with corresponding results of various computa- 
tions, we adjust the main stellar modeling parameters: mass, age 
and metallicity. In Fig.[3j the rectangular error boxes correspond 
to the values and accuracies of the T e ff and L parameters quoted 
in Table [5] The values of the radii reported in the present work 
select very narrow diagonal sub-areas in these error boxes. The 
new measurements of radii are thus particularly discriminating 
for the models, as previously noticed by Thevenin et al. (2005 ) 
and Creevey et al. ( 120071 ). 

We adopt an initial helium content of Y,„, =0.265 which cor- 
respond to a slightly metal poor disk star. We then tried to fit 
the evolutionary models within the error boxes of both stars 
with the masses proposed by Gorshanov et al. (2006| of 0.74 
and 0.46 M . In this process, we explored a range of abun- 
dance Z/X by changing the observed value by +0.15 dex and 
Y by +0.01 dex. We also tried varying the masses by ±0.05 M . 
No combination of the three parameters produced evolutionary 
models reaching the error boxes of the HR diagram, that are 
strongly constrained by the radius. The variation of the a param- 
eter of the convection depth did not help either. From this we 
conclude that the masses proposed by Gorshanov et al. ([2006 ) 
are not reproducible by our modeling within +0.05 M . For this 
reason, we decided to use as a starting value the older mass de- 
termination by Walker et al. ( [1995] ) of 0.67 and 0.59 M Q . The 
convergence of our models towards the observations is signifi- 
cantly better and we could refine these masses to the values listed 
in Tabled 

We selected as the most plausible models those satisfying 
first the luminosity and radius constraints and second the effec- 
tive temperature constraint. The corresponding parameters are 
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Fig. 3. Evolutionary tracks in the H-R diagram for 61 Cyg A 
(top) and B (bottom). The labels indicate the age in Gyr rela- 
tive to the ZAMS. The rectangular box represents the classical 
L-T e ff error box, and the diagonal lines represent the radius and 
its uncertainty. 

given in Table [5] The models of 61 Cyg A and B converge si- 
multaneously to the radii-limited uncertainty boxes for an age of 
6.0 +1.0 Gyr. This is significantly older than the 2.1-1.9 Gyr es- 
timate of Barnes (120071) . that was derived from the measured ro- 
tation period of the stars (gyrochronology). The chromospheric 
age quoted by this author of 2.4-3.8 Gyr also appears lower than 
our value. 



5. Asteroseismic frequency predictions 

Late-type bright binaries are rare but represent, once analysed 
with asteroseismic constraints, an excellent challenge for inter- 
nal structure models to derive the age and helium content in 
the solar neighborhood. With the improved determinations of 
the fundamental parameters of the 61 Cyg system, we propose 
predictions of pulsation frequency separations useful for astero- 
seismic diagnostics. To date, no detection of pulsation has been 
observed in these two stars. Our purpose here is not to predict 
which modes are excited but rather to compute a broad range 
of eigenfrequencies v„j with radial orders n > 10 and degrees 
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Fig. 4. Diagram representing the small separations Sv n> ( (vertical 
axis) as a function of the large ones Av„/ (horizontal axis) for 
the two components 61 Cyg A (A) and B (•). We consider here 
only (6v„fi, Av„ o) (solid line) and (<5v„,i, Av nj i) (dashed line). The 
successive radial nodes n > 10 are indicated on the curves. 

I = 0, 1, 2, 3 up to the cut-off frequency, so that they will cover 
possible future detections and eventually help to constrain fu- 
ture models of these stars. For a thourough review of the astero- 
seismic concepts discussed in the present Section, the interested 
reader is referred to the excellent review by Cunha et al. ( 120071 1. 

The frequencies v„j are calculated with a standard adiabatic 
code for stellar pulsations (e.g. Unno et al. 119891 1. Instead of 
looking at the absolute values of the frequencies, we consider 
the so-called small Sv„i = v„/ - v n -\^+2 and large separations 
Av„/ = v n j - v n -ij which are the relevant quantities to constrain 
stellar properties (e.g. Christensen-Dalsgaard 1984). The former 
gives information about the stellar evolutionary status of the star 
since low degree modes have their inner turning points close to 
or in the core of the star. It is particularly dependent on the gra- 
dient of the sound speed in the core which changes as the star 
evolves (e.g. Gough |T9861 l. The large frequency separation Avq 
roughly equals: 




with C j the sound speed. This quantity is a measure of the propa- 
gation time inside the star. It is proportional to the mean density 
which is then a strong constrain on the mass once the radius is 
determined by interferometry (Cunha et al. 2007). We plot these 
separations in Fig.|4]in a (<5v„/, Av„/) diagram for the two com- 
ponents 61 Cyg A & B. The small and large separations show 
the same behaviors with the radial nodes for the two stars. The 
shifts between Av n j for the two stars reflect the difference of 
their mean densities. It is interesting to note the change of sign 
of 5v„ ( for large values of n. 

6. Conclusion 

We presented high accuracy interferometric measurements of 
the angular diameters of the two nearby stars 61 Cyg A & B: 
<9 LD (A) = 1.775+0.013masandf? LD (B) = 1.58 1+0.022 mas, cor- 
responding to photospheric radii of R(A) = 0.665+0.005 R Q and 
R(B) = 0.595 ± 0.008 R Q . We computed CESAM2k models that 
reproduce these radii as well as the other observed properties of 



the two stars. 61 Cyg A & B appear as very promising targets for 
future asteroseismic studies, and we also derived asteroseismic 
frequencies potentially present in these two stars. The detection 
of oscillations would bring important constraints to stellar struc- 
ture models in the cool, low-mass part of the HR diagram, where 
convection plays a central role. However, in the absence of mea- 
sured asteroseismic frequencies, it appears difficult to go beyond 
the present modeling of the binary system. The main reason is 
that the mass is not constrained sufficiently well by the long pe- 
riod astrometric orbit. We encourage asteroseismic groups to in- 
clude these two stars in their observing programmes, as the mea- 
surement of seismic parameters (in particular the mean large fre- 
quency spacing) will bring a decisive constraint to the mass of 
these stars. We also found that the value of the mixing length 
parameters for both stars is very well constrained by the radius 
once the other parameters of the model are fixed, in particular 
the mass. We reach the same conclusions as Yildiz et al. (2006) 
that the a parameter is mass dependent at least for solar abun- 
dance stars. The future Gaia mission (Perryman 2005 ) will open 
the way to a complete calibration of several thousand binaries. 
Such a large sample will give us a completely renewed view of 
convection in stars (see e.g. Lebreton |2005l ). 
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